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The evolution a network of vortex loops due to the fusion and breakdown in 
the turbulent superfluid helium is studied. We perform investigation on the 
base of the "rate equation" for the distribution function n{l) of number of 
loops in space of their length I. There are two mechanisms for change of 
quantity n{l). Firstly, the function changes due to deterministic process of 
mutual friction, when the length grows or decreases depending on orienta- 
tion. Secondly, the change ofn{l) occurs due to random events when the loop 
crosses itself breaking down into two daughter or two loops collide merging 
into one larger loop. Accordingly the "rate equation" includes the "collision" 
term collecting random processes of fusion and breakdown and the deter- 
ministic term. Assuming, further, that processes of random colliding are 
fastest we are in position to study more slow processes related to determin- 
istic term. In this way we study the evolution of full length of vortex loops 
per unit volume-so called vortex line density C(t). It is shown this evolution 
to obey the famous Vinen equation. In conclusion we discuss properties of 
the Vinen equation from the point of view of the developed approach. PACS 
numbers: 6140.Vs 98.80.Cq 7.37.+q 



1. INTRODUCTION AND SCIENTIFIC BACKGROUND 

In spite of very long history the theory superfluid turbulence is very far 
from the more or less completeness ( 1 , 2 , 3 ). The reason to this is an incredible 
complexity of this problem. Indeed we have to deal with a set of objects 
(vortex loops) which do not have a fixed number of elements, they can split 
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and merge in processes of reconnection. Thus, some analog of the secondary 
quantization method is required with the difference that the objects (vortex 
loops) themselves possess an infinite number of degrees of freedom with 
very involved dynamics. Clearly this problem can hardly be resolved in the 
nearest future. As for direct numerical simulations, which remain the main 
source of information about this process 4 , 5 , 6 , 7 there are also many problems. 
It is clear that any progress can be achieved only if one can essentially reduce 
the number of degree of freedom of each of the loops. One way elaborated in 
context of lambda-transition ( 8 ) is to imagine vortex loops as a set of rings 
of different sizes and to take their radius as the only degrees of freedom. 
In paper of author (see 9 ) there is offered another way - to think of vortex 
loops as randomly walking chains. The physical ground for this supposition 
is the following one. During evolution the loops undergo huge amount of 
collisions with other loops (or with itself) with consequent reconnection. As 
a result each of the loops consists of many uncorrelated parts of mean size 
£o so called elementary step. A bit more precisely an average loop consists of 
smoothly (but randomly) connected arcs of mean radius £o- This fact allows 
to consider vortex loop as a random walk, and to describe its statistics with 
help of some generalized Wiener distribution. Parameters of this probability 
distribution functional can be regarded as degrees of freedom of vortex loop. 
This conception allowed to fulfil a series of studies on the vortex tangle 
properties, e. g. to calculate its momenta, energy and spectrum of energy. 
In the present work we demonstrate how this point conception allows to 
derive analytically the famous Vinen equation for evolution of the vortex 
line density C(t). 



2. KINETIC EQUATION 

Following the work 10 we introduce the distribution function n(l,t) of the 
density of a loop in the "space" of their lengths. It is defined as the number 
of loops (per unit volume) with lengths lying between I and I + dl. There 
are two mechanisms for change of n(l,t). One of them is the deterministic 
process. In superfluid helium quantity n(l, t) changes due to interaction of 
vortices with the normal component. This processes is more or less clear 
and can be modeled by relatively simple equations of motion. Other reasons 
for change of quantity n(l,t) are the random reconnection processes. We 
discriminate two types of the reconnection processes, namely the fusion of 
two loops into the larger single loop and the breakdown of single loop into 
two daughter loops. The kinetics of the vortex tangle is affected by the 
intensity of the introduced processes. The intensity of the first process is 
characterized by the rate A(l±,l2, 1) (number of events per unit time and per 
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unit volume) of collision of two loops with lengths l\ and h and forming the 
loop of length I = l\+l 2 . The intensity of the second process is characterized 
by the rate of self-intersection B(l,l±,l2) of the loop with length I into two 
daughter loops with lengths l\ and I2. In view of what has been exposed 
above we can directly write out the master "kinetic" equation for rate of 
change the function n(l,t). 

^dT + ^¥^Ft = I J MhM,l)n(h)n(l2)5(l-h-h)dhdl2 

- J J A{h, I, l 2 )n(l)n(h)6(h - h - l)dhdl 2 

- J J A(l 2 , 1, h)n{l)n{l2)5{h ~h~ l)dhdl 2 
- J J B(l u l 2 , l)n(l)5(l -h- h)dhdh 

+ J J B{lMM)n{h)5{h-l-l 2 )dhdh 

+ J J B(l, hM)n{li)5{h - I - h)dhdh . 

Coefficients A and B are calculated in 11 

A{h,l 2 , 1) = bmVihh, B(h, h, I) = bsVddoh)-^ 2 . 

Here b m and b s are some constants of order of unity. In paper 11 on 
vortex loops in superfluid turbulent Hell there was offered b rn pa 1/3, 
b s pa 0.0164772. Quantity V\ is the characteristic velocity of lines by or- 
der magnitude equal k/^q (k is the quantum of circulation). The Brownian 
random walk approach fails for scales near £q, therefore usually this value 
appears as a low cut-off. In paper 12 there was shown that the collision 
term alone leads to steady solution expressed by formula 

n{l) = Cr 5/2 . (2) 

Solution (2) is nonequilibrium solution with the I -independent flux P of the 
length density L = I nil) , which is the length accumulated in loops of length 
lying between I and I + dl. Term "flux" here means just the redistribution 
of length among the loops due to reconnections. Quantity P is equal 

P = 6. 2775C 2 b m K/£ - 2. 772 5C7) s k£ " 5/2 - (3) 

The master equation (1) is the base for study of superfluid turbulence in 
the frame of conception of randomly walking chains. The aim of our present 
paper is to demonstrate that the vortex line density, total length per unit 
volume defined as C(t) = J n(l)ldl evolves in accordance with the famous 
Vinen equation. 2 
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3. VINEN EQUATION 

To show it let us multiply the kinetic equation (1) by / and integrate 
over all sizes 

dt J dt J 31 dt v ; 

Here the quantity P ne t is the net "flux" of length in I— space, which is absolute 
value of flux P expressed by relation (3). This rule is imposed because flux P 
always carries away the vortex line density C from the system, and different 
signs refers to direction of the cascade. First we treat the deterministic term 
in equation (4). Due to enormous number of reconnection it is quite natural 
to suppose that "equilibrium" solution (2) is reached much faster than the 
slow deterministic change of function n(l). Then we can use this solution to 
evaluate deterministic term. We will calculate the rate of change of length 
of each loop on the base of the motion equation of the line in local approach 
(see e.g. 4 ) 

= (3s' x s" + as' x (V ns - /3s' x s") (5) 

Here s' and s" are the first and second derivatives from position of line s(£) 
with respect to label variable £, which coinsides here with the arclength. to 
calculate, v/(£) is velocity of the line element. To calculate 31/ dt, we use the 
relation for the rate of change of length 381 /dt for some arbitrary element 
with length 51 (see.g. 4 ). Assuming for a while that the label variable £ is 
not exactly the arclength, we have 51 = |s'| <5£. Then the following chain of 
relations takes place. 

351 3\s'\5i |s'| 3\s'\6t s' ds'5£ 



dt dt Is'l dt Is'l dt 



s'v^e 



On the last stage we return to |s'| = l.Diffirentiating (5) and multiplying by 
s' we have after little algebra 

^ = (a(s' x s")V ns - a/3(s' x s") 2 )^ (6) 

The next step is to average expression (6) over all possible configurations of 
vortex loops. We do it with use of the Gaussian model of the vortex tangle 
elaborated earlier by author. 9 In accordance with this model 

<S' X S") = J 1 ((s , x g „ )2) = ((s „ )2) = 1 (?) 

Quanity c 2 is the structure constant introduced by Schwarz 4 ). It follows 
from (6) and from (7) that 351 /dt oc I. Substituting (7) into (averaged equa- 
tion (6)) and then into (4) and integrating by part we get the contribution 
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into dC(t)/dt from the deterministic term. 

, h |V„,I ,, 1 , ( On(l,t) 2 2J,|V„,I r , ag 

Now we have to treat the "flux" term in equation (4). Unlike paper 12 where 
the constant C in relation (2) was obtained from condition P = 0, here we 
obtain C from the normalization condition C(t) = J n(l)ldl = C J l~ 3 / 2 dl = 
2C I £o (We recall that £ is a low cut-off of the wholle approach and integral 
diverges on low limit ). Furthermore we consider consequently the collision 
and reconnection events to put the system in equilibrium (with respect to 
solution (2)) state much faster than the slow deterministic processes. This 
implies that parameters £o and Ii, the so called structure constants of the 
vortex tangle, have a time to adjust to their equilibrium values. This as- 
sumption is widely adopted and it was confirmed in numerical simulations 
in. 13 In particular the mean radius of curvature (playing a role of elementary 
step in our approach) is related to the vortex line density as £g = 1/2 c| C 
(see 4 , 9 ). By use of the said above we rewrite expression for flux (3) in form 
P = CfkC 2 , where the temperature constant Cf is 

C F = (1.569 46 m - 2.772 6c|b s ). (9) 

We named this constant in honor of Feynman who was the first person to 
discuss evolution of vortex line density due to the reconnection processes. 
We would like to recall that he supposed decay of the vortex tangle due to 
the cascade like breakdown of vortex loops with further disappearance of 
them on very small scales. Relation (9) shows that there is possible the in- 
verse cascade, which corresponds to the cascade-like fusion of vortex loops. 
Unfortunately our approach has too approximate character to do any strong 
quantitative conclusion. It is clear, however, that for low temperatures, 
where the vortex tangle is more kinky, correspondingly C2 is large, the quan- 
tity Cf is negative. This corresponds to the direct cascade in region of 
very small loops. On the contrary for high temperature lines are smoother, 
C2 is small, and Cf is positive, which implies that there is inverse cascade 
with formation of large loops. If we for instance adopt values for b m and b s 
and use for c| values offered by Schwarz (see 4 ), then we get Cf ~ —0.01 
for the temperature 1.07 K and Cf ~ 0.44 K for the temperature 2.01 K. 
Collecting contribution into dC(t)/dt (4) from both the deterministic and 
collision processes, and taking into account that P net = \C F \ kC 2 we finally 
have 

^ = \ah \V ns \ d - \aPc 2 C 2 - \C F \ nC 2 . (10) 
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Thus, starting with kinetics of a network of vortex loops, we get the famous 
Vinen equation. 2 Let us discuss meaning of various terms entering this equa- 
tion. The first, generating term in the rhs of the Vinen equation describes 
the grows of the vortex tangle due to mutual friction. The second term also 
connected to mutual friction, however this term is responsible for decrease of 
the vortex line density. This point of view coinsides with ideas by Schwarz 4 
who obtained the deterministic contribution into dC(t)/dt using a bit dif- 
ferent approach. The third term in the RHS of (10) is related to random 
collisions of vortex loops. It describes decrease of the vortex line density due 
to the flux of length carrying away the length from the system. Depending 
on interplay between coefficients b m and b s and parameters ci the flux can 
be either positive or negative. Negative flux appears when break down of 
loops prevails and cascade-like process of generation of smaller and smaller 
loops forms. There exists a number of mechanisms of disappearance of rings 
on very small scales. It can be e.g. acoustic radiation, collapse of lines, 
Kelvin waves etc. Thus, in this case the situation is fully coincides with the 
scenario proposed by Feynman. 1 The case when the flux is positive is less 
clear. Positive flux implies the cascade-like process of generation of larger 
and larger loops. Unlike previous case of negative flux, there is no appar- 
ent mechanism for disappearance of very large loops. There are possible 
the following possibilities. The first one is that very large loops, whose 3-D 
size is comparable with size of the volume (note that length of these loop is 
much larger because of the random walk structure) can annihilate fully or 
partly on boundaries. Then the picture reminds the Feynman scenario but 
with the inverse flux. The second possibility is that parts of large loops is 
pinned on the walls and the whole approach representing the line to have 
the random walk structure fails. Finally it should not be ruled out, that a 
state with lines pinned on the walls and stretching from wall to wall is just 
a degenerate state of the vortex tangle. Some of numerical investigators 4 , 6 
report on this situation. 
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